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Abstract 

We consider a class of multi-matrix models with an action which is 0(Z?)-invariant, where D is 
j^H. the number of N x N Hermitian matrices X^, /i = 1, . . . , D. The action is a function of all the 

elementary symmetric functions of the matrix = Tr{X^X v )/N . We address the issue whether 
the O(D) symmetry is spontaneously broken when the size N of the matrices goes to infinity. 
The phase diagram in the space of the parameters of the model reveals the existence of a critical 
. boundary where the 0{D) symmetry is maximally broken. 

S3' 



OUTP-02-29P 



* vernizzi@thphys.ox.ac.uk 
' j.wheaterl@physics.ox.ac.uk 



1 Introduction 



Over the last twenty years several multi-matrix models have been considered for the description of 
a wide range of physical systems, from Statistical Physics to QCD or Quantum Gravity ||, ||[ |J. 
Although an analytic solution is generally not as easy to achieve as for the single-matrix models, a 
remarkable number of successes and results have been obtained so far @. A general feature of one- 
matrix models is that they possess an internal global symmetry under some gauge group (e.g. U(N) 
invariance, where N is the size of the matrix) which determines much of the universal behaviour in 
the large N limit. This global symmetry is present also in all the most relevant multi-matrix models 
(Ising model on random lattice [||, fjj, the Q-state Potts model [||, ||, 10, 11, 12, 12], chain of matrices 
models for coloring problem [20, 22, 23, 24], vertex models p5L | 



ie, 17 




the meander model |30, 31], the 0(n)-model and some generalizations of it |3_ 
IH U, and several others j4|, ||, |3| || ||, ||, ||, ||, The list is not complete). However, they 
do not usually have any further symmetry, except for the 0(n) model and its generalizations where the 
whole set of matrices transform as a 0(n) vector. The symmetry of these models is then U(N) x 0(n). 
Recently a new class of multi-matrix models have been introduced in the framework of Superstring 

1 and 



Theory and M-theory and the two main representative are the so-called IKKT-model |5lL 52 



the BFSS model 54]. They are proposed to be a non-perturbative definition of type IIB Superstring 
theory and M-Theory, respectively. In particular the IKKT-model is just one element of a bigger 
class of matrix models, called Super Yang-Mills Integrals (for an introduction see 
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The latter are characterized by carrying several (super) symmetries and they are obtained from the 
complete dimensional reduction of D-dimensional SU(N) Super Yang-Mills theories. These integrals 
also might provide an effective tool for the calculation of the bulk Witten index of a supersymmetric 
quantum mechanics theory [pT 



J |6_L],|62|, |6J. 

One consequence of having several symmetries is the existence of flat directions in the action of the 
model. They are potential sources of divergences when evaluating the integrals. The precise domain of 
existence of all the Yang-Mills Integrals with and without supersymmetry and for all the gauge groups, 

(after numerical and analytical studies for small gauge 
I [70|, [Q|, [72], |73|, [74], The existence 



has been rigorously determined in [56, ^57j, |5£ 

and for large gauge groups in 



groups m 



361, 67, 




of such flat directions affects not only the convergence properties of the Super Yang-Mills integrals, 
but also the behaviour of all the correlation functions and of the spectral density asymptotics. During 
the last few years it has been claimed that the "rotational" 0(D) symmetry (where D is the number 
of matrices) might be spontaneously broken in the large N limit |76]| . This issue has been analyzed 

and a possible 



in a series of analytical and numerical studies |6^, [7(], |7T], 72, 73, 74, 77, |78|, |8( 
mechanism for having such a spontaneous symmetry breaking has been proposed in 

The basic idea relies on the fact that these integrals contain fermionic degrees of freedom (i.e. 
matrices with Grassmannian entries) in such a way that the action is a complex number in general. 
However the action is a real number for lower-dimensional "degenerate" configurations (i.e. when 
the matrices are linearly dependent). Therefore, when summing over all possible configurations in 
the partition function the rapid oscillations of the complex action might enhance lower dimensional 
configurations in the large N limit. In order to shed light on this mechanism, a class of simplified 
fermionic multi-matrix models having a complex action (and the same 0(D) x U(N) symmetry) has 
been studied in [79|. In that case, the symmetry breaking actually occurs, and it is shown to be 
a consequence of the fact that the action is complex. Also, the results in j8^[ give indications of a 
spontaneous symmetry breaking in the IKKT-model. However, the actual mechanism for having such 
a behaviour (if confirmed) remains an open question. 

In this paper we address the question of whether a complex action is necessary if there is to be a 
spontaneous breaking of the 0(D) symmetry at large N. The action of the Super- Yang Mills Integrals 
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is complex in general but it also has flat directions. These two features have a quite different origin. 
The former is a consequence of the particular choice of the structure of the spinors together with the 
signature of the D-dimensional "space-time" in consideration. The latter arises because the action 
is made up of commutators or logarithms of fermionic determinants (or Pfaffians) (which are there 
ultimately as a consequence of having an highly symmetric theory). Since it happens that along the 
flat directions the action becomes real, it is not clear whether the spontaneous symmetry breaking 
is a consequence of the complexity of the action, or its flatness properties. A definite answer to 
this question would be given by a complete analytic solution of real-action models such as the Super 
Yang-Mills integral in four dimension, or its bosonic version (at any D) in which the fermions are 



suppressed. However only numerical simulations are available so far. The results of [82] suggest that 
there is no spontaneous symmetry breaking in the pure bosonic Yang-Mills integral. About the 4D 
Super Yang-Mills integral there has been some dispute [^, frl] whether there is symmetry breaking 
or not, and about which is the most reliable order parameter to use in that case (for a review see 



We decide then to focus our attention on building-up a multi-matrix model with a real positive 
semi-definite action made of standard Hermitian matrices ("bosonic"), but which allow a wide class of 
possible "degenerate" configurations. In this paper, we shall introduce a multi-matrix model sharing 
the same O(D) x U(N) symmetries, but with a real positive weights and without any Grassmannian 
degrees of freedom. This action allows many degenerate configurations and we will find that they can 
affect the symmetry of the model at large N. This fact is an indication that the exact mechanism 
which could be at the origin of a possible spontaneous symmetry breaking of rotational symmetries 
in Super Yang-Mills integrals deserves further studies. 

The paper is organized as follows: in Section |2| we define our multi-matrix model. It is based on 
all the elementary symmetric functions of the eigenvalues of the two-point correlation matrix, and it 
is manifestly O(D) invariant at finite N. The model contains a number of coupling constants which 
controls the role of the various elementary symmetric functions and the interaction among them. We 
study the behaviour of the model in the space of such parameters. In particular we solve the model in 
the simple and illuminating case where only two basic elementary symmetric functions are involved, 
i.e. the trace and the determinant. This case is simple enough for carrying explicit calculations at 
large ./V by means of a saddle-point method. In Section |3| we consider the more general case where all 
the elementary symmetric functions are present. There we show how the model is stable under such 
a generalization and that the 0(D) symmetry of the system holds everywhere except on a critical 
boundary where the symmetry is maximally broken. Finally, Section Q is devoted to our discussions 
and conclusions. For the sake of completeness, the appendix contains the calculation of a Jacobian we 
make use in Section ||. 

2 The model 

Let us consider a set of N x N Hermitian matrices {A^}, fj, = 1, . . . , D. The corresponding two-point 
correlation matrix 

T, u ^±Tr(X,X v ) (2.1) 

is a D x D real symmetric positive semi-definite matrix, with eigenvalues ti > • • • > tjj > 0. From 
the definition fl2.1| ) we see that if A M — > A^ = Q fJiU X u where Q G 0(D), then T transforms as 
T' = QTQ T . A straightforward consequence is that all the eigenvalues tn of T are 0(D) invariant 
quantities. Moreover, the matrices {A^} are linearly dependent iff some of the eigenvalues of the 
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correlation matrix T^ v are identically zeroQ. More precisely, a good indicator of the degree of non- 
degeneracy of the matrices {X^} is r(T) = rank(T), i.e. the number of non-zero eigenvalues of the 
matrix T. The most general action which is 0(D) invariant and is a function of the variables only, 
can be expressed in terms of the elementary symmetric functions Ck of the variables {t^}, k = 0, . . . , D. 
We recall here that the fc-th order elementary symmetric function of the variables {t^} is defined as 
the products of k distinct variables 

Ck = Vl^M2 • • • Vfc • (2-2) 

Ml<^2<---<Mfc 

(we omit the explicit independence of c^). It is well known that the Ck can be obtained from the 
expansion of the characteristic polynomial of the matrix T 

C(z) = det(I DxD + zT) = c D z D + c D - X z D - 1 + ... + <%. (2.3) 

All the Cfc are non-negative, as the matrix T is positive semi-definite. In particular one has 

1 ° ( 1 \ 

ci=trT=-^Tr(A^), c D = det I -Tr X,X u j , 

where we use the symbol "Tr" and "tr" to indicate the trace over N x N and D x D matrices, respec- 
tively. 



The partition function we consider in this paper is 



2\a\ 



D 
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-JvTtE£=x*2 



D 



a k N 2 



(2.4) 



k=l 



where ctk are real parameters. Eq. (2.4) is manifestly 0(D) invariant. This symmetry is not to be 
confused with the usual U(N) "internal" symmetry, which still holds for this model. In fact 2 [a] is 
invariant under X^ — * UX^U^ , for all fi, with U € U(N). The region of existence of this model as a 
function of the real parameters at will be determined later in this Section. Here we just emphasize 
that the argument of the matrix-integrals is always real and positive semi-definite. Moreover, another 
feature of eq. ( |2.4| ) is the existence of "flat directions" . They correspond to configurations where the 
matrices {X^} are linearly dependent, i.e. such that some of the symmetric functions Ck are identically 
zero. The convergence properties of the integral ( |2.4[ ) for large values of the entries of the matrices 
are mainly guaranteed by the presence of the Gaussian weight, but not completely. In fact, the flat 
directions contain non-integrable singularities (with some analogy to the case of the Yang-Mills inte- 
when some of the parameters are too negative. An exact bound in the space of 

is presented in eq. 



grals 
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|3.q). At finite N the average 



the parameters {ctk} for the existence of eq. (^ 
eigenvalues (tn) of the matrix T are all equal, because of the 0(D) invariance of eq. (2.4). However at 



large N this may no longer be the case, and our aim is to see whether the 0(D) rotational symmetry of 
the model can be spontaneously broken when N — > oo. In this context, we define also the dimension- 
ality d of a configuration of matrices {X^} as the number of non-vanishing eigenvalues of the average 
correlation matrix (T). Of course, at finite N one always has d = D. A possible way for probing 
O(D) symmetry breaking is to introduce an explicit symmetry breaking term before taking the large 

1 A short proof: if {X^} are linearly dependent, then 3??^ not all zero such that yj M rj^X^ = 0. Therefore T^r^v = 0, 
i.e. T has a zero eigenvalue. On the other hand, if ~Y^ v Tiivi) v = 0, then tr [(53 ^"m^m) 2 ] ~ Yin, v Vi^T^rt^ = which 
implies Y^u^^ 7 !^ = 0. 



4 



N limit. We do this by modifying the Gaussian weight in eq. (^4[) e~ N ^ ^ r ^) — > e ~ N T.^ A pTr ^ 
where the variables < Ai < A2 < • • • < Xd maximally break the O(D) symmetry of the model (in 
analogy to ||79||). After taking the large A limit, we shall remove the symmetry breaking term by 
taking the limit A^ — ► 1, V//. If (t^) — > for different directions /i then there is spontaneous symmetry 
breaking of the O(D) symmetry. 

We start with the simple case where Q2 = «3 = • • • = otD—i = 0- The partition function reads 

Z[a,A}= [ HdX^e-^^^' (txTr N2 (detT) aDN2 (2.5) 

It is convenient to introduce the matrix A = Suv^ui so that the partition function can be written as 
,Z[a, A] = /n^Li^A^ exp(— A 2 So) where the action is 

S [T,a,A) = tr(AT) - a x log tr T - a D log det T . (2.6) 

The action Sq depends on all the matrices only through the matrix T: it is therefore natural 
to change the integration measure from the multi-matrix variables {A^} to the single-matrix {T^ u }. 
When A 2 > D we have (see the Appendix [A| for details) 



Z[a,K]=C N , D dTe -^,A] + ^lo g ( d etT) ) c v 

Jt>o — T1 D r 1 N - fc + 1 



where the integral is over all the D x D real symmetric positive-definite matrices, the measure is 

N 2 — D — 1 

cff = Jl^>j, dT^ u and the Jacobian of the transformation is proportional to det(T) 2 . The 
partition function now has the proper form for the study of the large A limit by means of the saddle- 
point (Laplace) method for the asymptotic expansions of multidimensional integrals. According to 
this method, the main contribution to the integral comes from a small neighborhood of the critical 
points, i.e. global minimum points in this case, of the action (we drop 1/A 2 sub-leading terms) 

S[T,a,A] =trTA-ailogtrT-di)logdetT, (2.8) 

where ao = ao + 1/2. The minima of the function S can be at the boundary of the integration region 
or at the interior of it. In the latter case the necessary stationarity conditions for having a minimum 
are {saddle-point equations) 

° -S[T,a,A]=\ fl 6 flv -a 1 ^-a D (T- 1 ) Vft = 0, (2.9) 



± fl>U 



for all 1 < v < n < D. Note that multiplying eq. ( |2.9| ) by T and taking the trace gives tr (TA) = 
«i + Dao- Since tr (TA) > we have to look for solutions of eq. ( |2.9| ) in the region of the parameters 
plane {«i, ao} 

ai >-Da D . (2.10) 

The condition ( |2.10| ) is actually a bound on the domain of existence of the model at large A. In fact 
as we have already announced, the integral in eq. fl2.7|) exists only when the parameters at satisfy 
suitable constraints, and eq. ( |2.10| ) is one of them. Namely, the integrand function in eq. (2/7) does 



not have singularities in the integration region, except perhaps at the integration boundaries. At large 
values of the entries of T the integrand function is regular and integrable for any value of a^, being 
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bounded by the exponential factor. However, the behaviour close to the origin can give non-integrable 
singularities. This fact is evident when passing to the eigenvalues {t^\ of T = OtO T . It yields 



poo D p 

Z[a,A] ~ / TT di M |A(i)| / dO 

JO ^ JO(D) 



-N 2 tv OtO T A 




ai N 2 




a D N 



2 , JY — J>— 1 



(2.11) 



where A(t) is the Vandermonde determinant Yln<i/(tp~ Io(D) * s ^ e integral over DxD orthogonal 
matrices, with Haar measure dO, t is a diagonal matrix with diagonal elements ii,... ,to and "~" 
means "up to a (irrelevant) proportionality constant". From eq. ( |2.11| ) we see that, first, in order to 
have an integrable singularity at each of the (D — l)-dimensional boundary where only one ta = 0, it 
has to be 

N 2 - D - 1 



a D N z + 



> -1 



(2.12) 



At large N this condition simplifies to = + 1/2 > 0. Secondly, by rewriting the integral in 



eq. ( 2.11 ) from cartesian coordinates into multi-dimensional spherical coordinates, one has that the 
radial integration exists if and only if 



(D - 1) + D ( p -') + Ar% + 



N z — D — 1 



N a,D ] D > —1 



(2.13) 



Note that there are no contributions from the integral over the orthogonal group: in fact it is finite 
and regular in t, since it is an integral over a compact domain of an analytic function in its variables. 
At large the condition (2.13) is fulfilled by a\ + (^ + a£))D > which is precisely eq. (2.10). In 
summary, the region of existence of the model at large N is 



V = {{ak} ■ a\ + Dan > and old > 0} , 



(2.14) 



and it is depicted in figure ||. We point out that the model at large N is well-defined and finite also 
on the boundaries of V, i.e. B\ = {ap = 0, a.\ > 0} and £>o = {«i = —Dan, old > 0}. 




OLD 



- -Dol d 

Figure 1: Phase diagram of the model in eq. ( |2.5[) . The shaded region is where the partition function 
is divergent. The wiggle line is the region where the model is one-dimensional. In the remaining region 
the model maintains the full O(D) dimensionality. 
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If Q.B > then we see immediately that the global minima of S in eq. (2.8) cannot be on the 
boundary of the integration region. Otherwise the matrix T would have at least a zero eigenvalue, 
that is detT = 0, and there eq. (|2.8| ) gives S — > +oo. Therefore, in this case the critical points must 
be in the interior of the integration region. Let us then solve eq. (2.9) for qd > 0. It is straightforward 



to see that any matrix T which is a solution of eq. fl2.9| ) has to be diagonal. Defining T = 8^ u t^ eq. 



T) reads 



an 



A, 



trT 



l,...D. 



This system of algebraic equations can be solved easily. First, by summing eq. 
an equation for x = trT /ax, 



(2.15) 

(2.15) over fi we get 



7 



old 



1 



D 

y — 



1 



trT 



(2.16) 



For any given real 7 and A eq. ( [2.16|) is a rational algebraic equation with D solutions in the variable 
x. All the solutions are real. In fact, by writing the real and imaginary part of x = x' + ix" and using 
the fact that 7, A 6 R, it yields x" = 0. Among such D real solutions, we have to pick up the ones 
that make t^> because T has to be a positive semi-definite matrix. From eq. ( |2.15 ) we obtain that 



trT 



1 

Ai 



> ^- for ot\ > 



x = < 1 for ai < 



(2.17) 



which is satisfied by only one solution in each case. Namely, for a\ > is 7 > and the solution is 
the largest possible one (the one greater than 1/Ai) whereas for a\ < is 7 < and the solution is 
the one with x < (see figure ||) . 



7 > 




7 < 




Figure 2: Graphical representation of the r.h.s. of eq. ( |2.16| ) as a function of x = The circles 
indicate all the solutions. The black circles indicate the acceptable solutions. 



For any given point {ai,an > 0} in the interior of the parameter space T>, if we remove the 
symmetry breaking terms by taking — > 1 V/x, then the (unique) solution of eq. ( 2.16j ) is x = 1+D/j, 
i.e. trT = a\ + ar>D. Inserting this value in eq. (2.15) we read that in the large N limit all the 
eigenvalues are equal to (t^) = ot\/D + old- Hence, we conclude that in the region inside V with 
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old 7^ the model has a phase with dimensionality d = D and the 0{D) symmetry is preserved, as 
expected. In such a phase, the free energy T = — -™ logZ reads (with all A M = 1) 

T = (ai + Da D ) [1 - log (ai + Da D )\ + Da D log D . (2.18) 



On the boundary £>o of 2? where <X\ = —Dold the solution of eqs. (2.15) and ( |2.16 ) gives T = 



However, it would be wrong to conclude that the dimensionality is d = there, because actually in the 
limit old — > —ol\/D one has T — > with d — > -D. That means that there is no spontaneous symmetry 
breaking on the boundary i5o- 

Let us consider now the final case of the boundary B\ where old — > + first and A^ —* 1 for all 
/i afterwards. In such a limit one must consider a\ > in order to stay within the region T>, eq. 
fl2.14|) , and therefore 7 — ► +00 .0 According to eq. ( |2.16| ) and figure [l|, this fact may occur only when 
x — > 1/Ai. From eq. ( |2.15 ) we have 



OL~\ Oil 

lim trT=— - and ^ = -±5^. (2.19) 

a D ^0+ Ai Ai 

In other words, only one eigenvalue of the matrix T is not zero in this limit. Removing the symmetry 
breaking term by setting Ai — > 1 leads to a dimensionality d = 1, actually (i^) = aic^-a. That 
concludes our proof that the model in eq. (|2.5| ) has a maximal spontaneous symmetry breaking of 
O(D) symmetry whenever art — * + . 

It is interesting to notice that we could have had considered directly the case old = (and not just 
the limit old — * + ), because there the model at large N is well-defined. In fact, let old = from the 
very beginning in eq. (|2.8|). Then, as the A^'s are all different each other, S cannot have any minima 
in the interior of the integration domain (in other words, eq. (|2.9| ) do not admit any solution). Hence, 
the global minimum must be on the boundaries of the integration region, where some tn is equal to 
zero. Analyzing by inspections all the hyperplanes which constitute the integration boundary, one 
finds that the global minimum is a point on the line ti =t% = ... = to,ti > and it is precisely at 
h = olx/X\. Substituting this value in eq. fl2.8|) gives the free energy for the phase B\ 

T = ol 1 (\-\°^). (2.20) 

This expression (for A = 1) matches continuously with the free energy in the un-broken phase, eq. 
( p. 18 ) for old — > + . By taking derivatives of the free energies with respect to A^ we can compute the 



correlation functions, in particular the average of the eigenvalues, and the susceptibility 

d 2 T 



. . dT 



d\ad\ u 



= -N z (t ll U) conn . (2.21) 

A=l 



A=l "-"V^ 
In the broken phase B\, we get from eq. (jOCp and (gjg ) 

(t,j.) = aiS^i , Xiiv = -aiSiiiSxu ■ (2.22) 



which is of course consistent with eq. ( 2.19 ). The computation of the same quantities in the un- 
broken phase requires the knowledge of an expression of the free energy as a function of A^ (i.e. eq. 



( p. 18 ) is not useful for that). A general analytic expression seems not so easy to get since it needs 



the analytic solutions of the algebraic equation ( p.!6|) in a closed form, which is known to be an 
impossible task when the degree of the equation is large. However, we can proceed as follows. We 



2 In principle it would be possible to take the same limit with Qi < but then one necessarily would end up in the 
origin of the coordinates ai = qd = where the system is purely Gaussian. 
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already know the pattern of symmetry breaking from eq. ( |2.19| ), Hence we can restrict to the case 
where Ai < A2 = • • • = Xd without loosing in generality. In this case, eq. ( |2.16| ) is a second order 
algebraic equation which can be solved explicitly. We obtain then the free energy, its first and second 
derivatives w.r.t. A^ and in the limit A^ — > 1 they are 

. ai 1 (a\ _ \ f —aUD - 1) - a D D 2 if a = u ,..„, 

Note that the susceptibility is divergent as ~ 1/cen when — > . The singular behaviour of the 
susceptibility is again a signal of a criticality at ar> = 0, where the rotational symmetry is actually 
maximally broken down to one dimension. 



3 Generalization 



Let us consider now the more general case eq. (2.4) where all the symmetric functions are allowed 
(and not only c\ and cd, i.e. the trace and the determinant, respectively). Again introducing the 
symmetry breaking term = X^S^, < X\ < A2 < • • • < Xd, and following the same path of 
reasoning as in the previous paragraph, we have 

Z[a,A] =C N , D f d Te- N2s ^ T > a ^ , (3.1) 
Jt>o 

where the generalized action at finite iV is now: 

D D + l 

S N [T, a, A] = tr TA - ^ a k log c k + log c D , (3.2) 

k=l 

and 6ik = ak + \$k D- Let us first determine the region of the parameter space where the partition 
function eq. (3.1) exists. To that aim it is worthwhile to pass to the eigenvalues t\,t2, ■ ■ ■ , £d of T 
in the integral ( |3.1| ), as we did in eq. ( |2.11 ), thus obtaining a D-dimensional integral. The condition 
which prevents there being a singularity at the point where all the t^s are zero is 



DID -I) p/A, D(D + l)\ 
(D-l)+ { 2 +^ 2 (E fc ^- j >_1 (3 - 3) 

as one can see by passing to high-dimensional polar coordinates^. More generally, the integrand 
function does not have singularities on the p-dimensional hyperplanes where D — p of the variables 
tu's are zero if and only if 

(D - p - 1)+ (°-ri(^-l> +y Jf (*-,),>- <°-g(f + 1 A >-l, (3.4) 

\fc=P+l / 

for < p < D — 1. In the large TV limit, the conditions in eq. ( |3.4| ) relax to 

D 

Y,( k -P)ak>0, p = 0,...,D-l. (3.5) 
k=p+i 



3 The first term of eq. ( |3.3| ) is the contribution from the radial part of the polar measure, the second is from the 
Vandermonde, and the remaining terms are from the action. The integral over the orthogonal group does not generate 
any singularity. 
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In particular note that old > 0. We call T> the region in the parameter space {cifc} which is determined 
by the conditions in eq. (|3.5|), and from now on we shall consider only values of the parameters {a k } 
which belong to T>. Obviously, this is a natural generalization of the analogous region obtained in eq. 



The generalized action eq. ([3^) at large N reads 



D 



S[T, a, A] = tr TA - ^ a k log c k , (3.6) 



k=l 

and in the same limit the main contribution to the partition function ( |3.1|) comes from the global 
minima of S. Such minima can be in the interior of the integration region or on the boundaries of it. 
In the former case, the saddle-point equations are 

d r , _ 1 dc k 



S[T 1 a,A] = X^-Y j a k - 7 ^=0. (3.7) 

Any matrix T which is a solution of eq. ( |3.7| ) must be diagonal. In fact, taking the commutator of 
eq. ( |3.7l ) with T yields [A, T] = because T commutes with any other function of T. Writing the 
commutator in components reads (A M — X U )T^ U = 0, i.e. T is diagonal. Thus, letting T = S^ w t^, the 
saddle-point equations are equivalent to the following system of non-linear algebraic equations 

A At = ^a fc — fi=l,...,D. (3.8) 
k=i Ck m » 

The case where the absolute minima of the action S are instead on the boundary of the integration 
region can occur only if some parameters a k are identically zero. In fact, if all the parameters a k are 
different from zero, then the action is positively divergent when at least one tn is zero, and thus there 
cannot be any minima on the boundary. 

For the moment, let us restrict the discussion to the case where all the parameters a k are strictly 
positive. We call T> + C T> such a region of the parameter space. It is straightforward then to show that 
in V + the system in eq. ( p.8[ ) has only one real positive solution (that is a set of {t\ > 0, . . . to > 0} 
which fulfills eq. ( |3.8|) ), and it is actually the single global minimum of eq. ( |3.6[ ). In fact, the linear 
combination tr TA = X^t^ and all the elementary symmetric functions c k are multilinear (/c-affme) 
functions in the variables ii, . . . jio, as one can directly see from the definition ( |2.2| ). As such they 
are convex functions. Also the function — log(a;) is convex for x > 0, and therefore the action S in 



eq. (3^) is a convex function, being a finite linear combination with positive coefficients of convex 
functions. Moreover, we show that S is also bounded from below. In fact, we can prove it by using 
the following inequality 



D D 



^ a k log c k < ^2 kdl k log ci . (3.9) 



k=l k=l 



The proof of the inequality (3J)) is by induction. For D = 1 it is an identity. Let us suppose that eq. 



( |3.g| ) is valid for D — 1. Therefore we have 

D D-l 



^ a k log c fc < ^2 kdt k log ci + a D log 

k=l k=l 

D-l / \ D D 



CD CD-I Ci 



CD-I CD-2 C 



u—s. / \D ±J 

< ^2 ka k log ci + a D log ( — j =^fca fe logci, (3.10) 
k=i VC °^ k=\ 
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where we used repeatedly Newton's inequalities c\ > Ck-iCk+i for 1 < k < D — 1, in the form 
Ck+i/ck < C/t/cfc-1, and the fact that old is positive in V. By applying the inequality ( |3.9| ) to the 
effective action eq. (|3.6D we get 



D 



S[T,a,A] > Aid-^fcajfclogci, (3.11) 



fc=i 



because tr TA > Ai X^=i Since the function ax — blog(x) > 6(1 — log |) for any o, 6, x real and 
positive, we finally obtain a lower bound for the action 

S[T, a, A] >a(i — log , (3.12) 



where A = Y^k=i ^ k * s positive in T>, as follows from eq. ( |3.5| ) with p = 0.0 

All the above shows that when {a^} £ D + the action S is continuous, lower bounded and convex 
in the integration region. From the additional observation that the action is linearly divergent when 
any is large and logarithmically divergent when any £ M is close to zero we conclude that necessarily 
the action has one and only one global minimum, and it must be in the region > 0, V/i. We call 
such a minimum i = {ii, . . . , to}, > 0. 

The large N limit of the model is controlled by the behaviour of i as a function of d&. In the 
following we enumerate a series of properties of t. To that aim is worthwhile to recall two useful 
properties of the elementary symmetric functions |9l], |92j. First, the fe-th order symmetric function 
Ck can always be decomposed as the sum of a i^-dependent part and a ^-independent part: 

c k = Vi-i + <#° , (3-13) 

where we defined = Cfc| t =0 , i-e. the k-th elementary symmetric function of {ii, t2, ■ ■ ■ , tp} omitting 
t^. Note that d^c^ = c£-v Second, the following equality holds 

D 

=Ck+i, k = 0,...,D-l. (3.14) 

Let us see now what consequences these properties have on t. 
1. The solution t of eq. (|3.8[) is upper bounded by 

D r U) D 



= ^ a k M k 1 < ^ a k , v> = 1, 



fc=i ft fc=i 
because from eq. ( 3.13 ) > i^c^ 1 . 
2. The solution f of eq. (|3.8[ ) is lower bounded by 

i D F M 

^t^*"—-^ (3 - 15) 

^ fe=l K f 1 

because all the terms in the sum are non negative and i fl cf D _ 1 / 'cd = 1. Therefore old has to 
go to zero for tn — * 0. Note that this condition means that when old > there cannot be 
any spontaneous symmetry breaking at all, since none of the eigenvalues is vanishing. In other 
words, if there is a phase transition, it must be on the plane old = 0. 



4 Note that the lower bound in eq. ( |3.12| ) is actually valid everywhere in T>, and not only for 6 T> + as our proof 

does not rely on such a restrictive hypothesis. 
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The minima are in general monotonic with respect to fi. Subtracting two equations of the 
system (|3.8[) gives 



D 



A u — A;. 



k=2 



-fc-2 



(3.16) 



and then the ordering Ai < A2 < . . . < \r> implies t\ > ti > . . . > if). On the other hand , from 
eq. ( |3.16| ) follows also £, = i u if and only if A^ = X v i.e. when the symmetry breaking terms 
are removed. We deduce that at any point of the region T> + , the dimensionality of the system 
is d = D and the original 0(D) symmetry is fully preserved. In this case we obtain (with all 

A„ = l) 



<*/*> = 4> T = A(l-logA) 



D 



D 

£ 

k=l 



ak log 







D k ( 





(3.17) 



with A = Y2k=l k a k 



4. Let us consider now the limit ocd — ► 0, while keeping all the other a.k<D fixed. In this limit the 
free energy has to be continuous either there is a symmetry breaking or is not. Its limiting value 



D-l 



•F d =A'[l- log A') log 



fe=i 



, i.e. 








w ( 


:?)] 



(3.18) 



with A = *^2k = i kak- If there is symmetry breaking then tp 



(it is the smallest eigenvalue) 

but no other i„ can go to zero. This is because if there are at least two iD-i,tr> — > then 
cd-i — > and eq. ( |3.8| ) would be inconsistent in the limit (l.h.s. is finite whereas r.h.s. is 
infinite). The free energy for a (D — l)-dimensional broken phase (with old = 0) would be: 



D-l 



Td-i =A'(1- log A) -J2 & * lo £ 



fc=l 



1 



(D - If 



D - 1 

k 



(3.19) 



In general Td < Fd-i with the equality only for D = 2, or D > 2 and S2 = ■ ■ • = «D-i = 0. 
We conclude that there is not spontaneous symmetry breaking when old — > 0, unless for D = 2, 
or D > 2 and 0.2 = ■ ■ ■ = ocd—i = (which is actually the case we considered in Section Q). 

It remains to consider the "wedge" region D/T> + of the phase space, where some of the ak are 
negative. In this case the action S is no longer a convex function, but it still possible to prove that 
it has only one global minimum. The proof goes as follows. First of all, S is still lower bounded by 
the same bound as in eq. ( 3.12j ), and it is divergent towards +00 at the boundaries of the integration 
region, hence it must have at least one local minimum. Secondly, if S has more than one local 
minimum then the system of equations ( |3.8[ ) would have multiple solutions for a set of values of the 
parameters {a^}- We know already that when {a^} is in T> + the solution is unique, therefore there 
must exists a value {a' k } of the parameters where multiple solutions merge together into the unique 
one. This implies that the Jacobian det di(a)/da has to be singular or zero for that particular value 
of a'. However we show now that this is not possible. In fact, let us write the system of equations 
Q3,8| ) in the more compact form: 

A = G[t(a)] ■ a (3.20) 

where A = (Ai, . . . , Xd), G>fc[i] = c il\/ c fc an d a = . . . , old)- Equation ( |3.20D implicitly defines 
the vector function t(a) as a function of a. We take the total derivative of the component \x of eq. 
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( |3,20 ) w.r.t. a.{ and compute the determinant w.r.t the indexes fx, i of the obtained expression. One 

det {G[t{ f • 5) " det %g> = (-!)» det G[t(a)] . (3.21) 

The first determinant on the l.h.s. of eq. ( |3,21 ) is regular and not zero at {a'}, otherwise the eq. (|3.8[) 
would not admit any implicit solution but we know it must exist (because the existence of a global 
minimum). The determinant on the r.h.s. is detG = A(i)/ IljLi c i' where A(i) is the Vandermonde 
determinant. This expression is finite and it zero only if at least two eigenvalues t^,t u are equal each 
other and then, by means of eq. ( |3.8| ), it must be A M = A„ which is not possible by hypothesis. This 
ends the proof that for any given in the "wedge" regions the action has only one local minimum 
in the interior of T>, which is then also a global one. The qualitative behaviour of this critical point 
as a function of the parameters a goes as for the case in T> + . After removing the symmetry breaking 
terms A — ► 1, the critical point becomes completely symmetric in its variables and it corresponds to 
an un-broken phase with 0(D) symmetry. 



4 Discussion and conclusions 

In this paper we have introduced a multi-matrix model where the Hermitian matrices X„ are inter- 
acting through all the elementary symmetric functions of the correlation matrix T^ u = Tr (X^X^/N . 
The main reason for the choice of such a model relies in its interesting features: first, it is manifestly 
0(D) x SU(N) invariant, and it allows the study of the issue of the spontaneous symmetry breaking of 
0(D) symmetry in the large- N limit. Secondly, the action of the model is real, positive definite and it 
does not contain any Grassmann variables. This is most useful for understanding what we can actually 
expect from a model without a complex action or rapidly fluctuating potentials. Understanding the 
effect of a complex action, which is a notorious difficult problem, requires also realizing first what 
could happen when it is not there. Third, it allows a number of possible "degenerate configurations" 
in the matrix integration measure and our aim is to understand their role in a scenario of spontaneous 
symmetry breaking. Finally, the model is considerably simple and can be solved analytically, being 
the interaction among the matrices only through the 0(D) "spatial" symmetry and not through the 
SU(iV) "internal" symmetry (for which there is just a Gaussian weight). We introduced a number 
of parameters which allows to tune the relative weight of the elementary symmetric functions of the 
model, and then we focused our attention on the phases of the model in the space of the parameters 
when N is large. This has been done in two steps: first in Section ^ by studying in full detail a 
simple case where only two symmetric functions are "switched on" (the trace and the determinant), 
and afterwards in Section || by considering the more general case where all the symmetric functions 
are present at the same time. In both cases we found that the 0(D) symmetry is broken only in the 
limit ao — > —1/2 for D = 2 or for D > 2 and oti = . . . = oid~i = 0. In these cases the dimensionality 
of the model collapses down to one dimension. 

The qualitative explanation of such a behaviour is simple. Degenerate configurations of the ma- 
trices such that the correlation function T^ v has zero eigenvalues, dominate the matrix integration in 
the large N limit, when the parameters of the model are tuned to a critical value. In particular the 
parameter old (which is coupled to the determinant, i.e. the elementary symmetric function most sen- 
sitive to "degenerate" configurations) is to be tuned to the critical value a£> = —1/2 for compensating 
an analogous "centrifugal" term coming from the Jacobian (see Appendix). At that precise value 
of otp = —1/2, the measure collapse down to one dimensional configurations, quite independently 
from the presence of other symmetric functions but the trace. This is most evident from the explicit 
solutions in Sections |2[ 
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The symmetry breaking mechanism of the model in this paper is therefore due to the existence 
of directions in the matrix integral along which the measure is identically zero. These directions 
are where the matrices are linearly dependent, with different degree of degeneracy. We learned also 
that the reality of the action does not seem to stop a generic Hermitian multi-matrix models with 
O(D) x SU(N) symmetry from having a spontaneous symmetry breaking of 0(D) symmetry when 
N is large. Of course this does not prevent other real-action multi-matrix models having different 
patterns of spontaneous symmetry breaking, nor does not say anything about the role played by a 
possible complex term in the action. For all these reasons our findings do not contradict the analysis 
of |6^, [7(], [n], [72], 73, 74, 77, 8^, 82]. It would be interesting to carry out the analysis contained in this 



paper to an extended version of the model where the coupling constants a are allowed to be complex 
numbers. The action would be complex then, and a different pattern of symmetry breaking seems to 
be possible. 

There are extensions of the model where the matrices are not Hermitian but real symmetric or sym- 
plectic. The only changes are in slightly different factors in the Jacobian (see Appendix) and they 
do not affect the large N results of this paper which still would hold in those generalized cases. We 
conclude by observing that the reason why we can solve this multi-matrix model is that the interaction 
among the matrices is only through the correlation matrix T^ u . For the rest the matrices are actually 
not interacting with the full internal SU(N) symmetry group, the interaction being just a Gaussian 
factor. In fact adding a quartic or higher order term to the action (i.e. terms like Tr X^X U X^X U and 
Tr A^A^) would probably change drastically this scenario, but it would also be more difficult to solve, 
as happens for multi- matrix models like the Yang-Mills integrals. 
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A Appendix 

For the sake of readability, in this Appendix we compute the Jacobian of the transformation in eq. ( |2.7| ) . 



It is a well-known result which has appeared several times in the literature, e.g. pa, pq, 87, pq, 



The technique we use here is similar to the one in [85]. The integral in eq. ( |2.4| ) is of the form 

D 

ndA^/(T[X]) (A.l) 

where / is a real function and the [/(A)-invariant integration measure for each Hermitian matrix is 
dX = Yl i=1 dXa Yii>j dReXij dlmXij as usual. First, by inserting the definition eq. ( EO| ) of the matrix 



T in the formula ( |A.1| ) by means of Dirac <5-functions, we can equivalently write 

D D 

1= dT f{T)J(T) , J(T) = / TT dX„ TT 5{T aP - -Tr X a X p ) . 
Jt >° J £=i «>p N 

The Jacobian J(T) can be evaluated by using the integral representation of the (5-function 

D D 



m = /n^n / 



dn 



H=l a>/3' 



2tt 
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where in the last equation we performed the Gaussian integral over the matrices X„, and we collected 
the elements £l^ u into a real symmetric matrix f2 (giving an additional factor from the measure). The 
real-symmetric matrix T can be diagonalized by an orthogonal matrix O, i.e. T = OtO T where t is 
a diagonal matrix, with diagonal elements > 0. We change the matrix variables — > W where 
W = O t CIO and we have dW = dfi, det(zW) = det(iO) and 



J(T) = C N)D f dW 



Jtr Wt 



det(iW) N2 / 2 



& itvW 



C „. D (detT)^ / iW^^ (A.3) 



where in the last equation we apply the transformation W^„ — ► W^ u / ^t^t u . The remaining T- 

independent integral is completely factorized and it is equal to 2 tt * / llfc=i -^C 2~~/ - Finally 

we obtain 

N ml n (2N 2_ D+1] 

J(T) = n^-D D („* (d6t T) ■ ( A - 4) 

The results in this appendix are valid for iV 2 > D, which is fine for the large N analysis of this paper. 
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